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Information that is stored in quantum-mechanical systems can be easily lost because of the inter¬ 
action with the environment in a process known as decoherence. Possible physical implementations 
of many processes in quantum information theory involve systems of identical particles, whence com¬ 
prehension of the dynamics of entanglement induced by decoherence processes in identical-particle 
open systems becomes relevant. Here we study the effects, and concomitant entanglement evolution, 
arising from the interaction between a system of two identical fermions and the environment, for 
two paradigmatic quantum channels. New entanglement measures are introduced to quantify the 
entanglement between the different parties, and a study of the dynamics of entanglement for some 
particular examples is carried out. Our analysis, which includes also the evolution of an entangle¬ 
ment indicator based on an entropic criteria, offers new insights into the dynamics of entanglement 
in open systems of identical particles, involving the emergence of multipartite genuine entangle¬ 
ment. The results improve our understanding of the phenomenon of decoherence and will provide 
new strategies to control it. 


I. INTRODUCTION 

Decoherence is a quantum phenomenon that plays an 
important role in connection both with the foundations 
of quantum physics and with its technological applica¬ 
tions m- Decoherence is also intimately connected to 
another key ingredient of the quantum world: quantum 
entanglement j3J. Indeed, the various effects associated 
with the phenomenon of decoherence are due to the emer¬ 
gence of entanglement between the system under study 
and its environment. The interaction between an imper¬ 
fectly isolated system and its surroundings, leads to the 
gradual disappearance of several quantum features ex¬ 
hibited by the system. These effects are at the core of 
the nowadays orthodox, decoherence-based explanation 
of the emergence of the classical world from quantum 
physics [suai 

Possible physical implementations of most of the pro¬ 
cesses in quantum information involve systems composed 
of identical particles 7j. Just to mention a few, sys¬ 
tems such as semiconductor quantum dots, in which 
charge carriers are confined in all three dimensions [8], 
or neutral atoms in magnetic [9] or optical unj micro¬ 
traps. However, contrary to what occurs in multipar¬ 
tite systems of distinguishable particles -in which much 
attention has been paid to analyse the dynamics of en¬ 
tanglement (mainly in qubit systems)-, the evolution of 
entanglement and decoherence in systems composed of 
indistinguishable particles remains a largely unexplored 
field. On one hand, the concept of entanglement in these 
kind of systems exhibits some differences from the corre¬ 
sponding concept as applied to systems consisting of dis¬ 
tinguishable subsystems, being perhaps more controver¬ 
sial. However, there is general consensus that in systems 
of identical fermions the minimum quantum correlations 


between the particles that are required by the indistin- 
guishability and the anti-symmetry of the fermionic state 
do not contribute to the state’s entanglement 0IIIM]. 
On the other hand, a resource theory based on local oper¬ 
ations and classical communication (LOCCs) is no longer 
suitable: Due to the necessary (anti)symmetrization of 
identical particles’ states, the notion of local operations 
does not apply anymore. 

The aim of the present contribution is thus to advance 
in the investigation of the dynamics of entanglement and 
decoherence in an open system composed of identical par¬ 
ticles. Specifically, we analyse the effects, and concomi¬ 
tant entanglement evolution, arising from the interaction 
between a pair of identical fermions and the environment. 

The article is organised as follows. Exchange- 
symmetry preserving transformations are discussed in 
section II, where it is evinced that the notion of local op¬ 
erations is foreign to identical-particle systems. Section 
III contains the preliminaries for the subsequent analy¬ 
sis of the entanglement distribution. First, we present a 
brief review of the definition and quantification of the en¬ 
tanglement between two identical fermions. Then we de¬ 
fine appropriate measures for the entanglement between 
one fermion and the rest of the system (fermion plus en¬ 
vironment), and also between one of the fermions and 
the environment. In section IV we investigate the dy¬ 
namics of the entanglement in an open system consist¬ 
ing of two identical fermions. We do this by considering 
two decoherence processes that are paradigmatic in quan¬ 
tum information theory. Our results open the possibility 
to study the efficiency of some previously developed en¬ 
tropic entanglement criteria, a task that is carried out in 
section V. Finally, some conclusions are drawn in section 
VI. 
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II. EXCHANGE-SYMMETRY PRESERVING 
TRANSFORMATIONS 

Consider a bipartite system S composed of a pair of 
subsystems a and b, immersed in an environment E. Ini¬ 
tially S is in an arbitrary state p a &(0), whereas E is as¬ 
sumed to be in a pure state denoted as |0) E . The initial 
density matrix of the complete system is thus given by 

p(0) = Pab( 0) ® |0> (0\ E . (1) 

At t = 0, S starts to interact with the environment E by 
means of a unitary transformation U = exp (—iHt/h), 
with H the total (S plus E) Hamiltonian. As a result, 
the effective evolution of the bipartite system ab writes 
as 

Pab(t) = Tr E (Up(0)U t) = ]T K,(t)p ab (0)Kl(t), (2) 

where the I\ M = (p\U\0) E are the Kraus operators 
associated to the transformation U, and {|/x) E } is an 
orthonormal basis of the Hilbert space He- Since 
TV p a b{t) = 1, the Kraus operators satisfy the relation 

Em K v K l = L 

We now assume that both subsystems a and b pos¬ 
sess the same dimension, say n. The n 2 -dimensional 
Hilbert space Hs can thus be decomposed into two sub¬ 
spaces, namely with diniH_ = n{n — l)/2, and H+ 
with diniH + = n(n + l)/2, that are spanned by basis 
vectors that are antisymmetric and symmetric, respec¬ 
tively, under the exchange of the subsystems a and b. 
Let {|'0fc ) s } be an orthonormal basis of the subspace 
n ±) where k= 1, 2, ...n(n± l)/2. If the initial state of the 
composite system S has a well-defined symmetry under 
the exchange a <-»• b, it decomposes as 

Pab(0) = J2 Pk U’k)(' , l J k\ab- ( 3 ) 

k 

Clearly, a necessary and sufficient condition for the 
evolved state (|2|) to preserve the initial symmetry is that 

[H,P ab ®I E } = 0, (4) 

with P a i, the subsystem-exchange operator. 

Up to now we have referred to E as the environment in 
which the bipartite system ab is immersed, without mak¬ 
ing further assumptions about its nature. In particular, 
E may also be a bipartite system, composed of indepen¬ 
dent environments A and B. Such decomposition of E 
allows for a local evolution, in which each subsystem a 
and b can couple independently (or locally) with its own 
environment A and B , respectively. A local evolution 
corresponds thus to a Hamiltonian of the form 

H = H aA + H' hn (5) 

where H aA = H[ ree + H f / ee + and similarly 

for H bB . In this case the unitary evolution factorizes 


as U=UaAU' bB , and the Kraus operators decompose as 
^=(>a,mb) = (am| U aA |(U> <S> {pb\ U bB |0 B ). If, on the 
contrary, E does not decompose into two independent 
environments but stands for the common environment 
of both (noninteracting) subsystems a and b, the total 
Hamiltonian writes as 

H = H[ ree + H' b free + H f E ree + H™* + . (6) 

This case corresponds to a global evolution, in which U 
cannot in general be decomposed as U aE U bE , and the 
subsystems a and b do not evolve independently. 

Now, it is easy to see that the condition 0 holds iff 
the evolution is global and symmetric, i.e., if H has the 
estructure ([h]) with H=H' for all the free and interaction 
Hamiltonians. Indeed, for 

H = Ht™ + H[ ree + H f E ree + H™ e + Htf, (7) 

it is immediate to verify that ( P a & (g> I E )H(P ab (g> 1 E )=H. 
The converse, that Q implies a global, symmetric evo¬ 
lution, can be verified assuming that the evolution is lo¬ 
cal and symmetric and showing that the corresponding 
Hamiltonian, namely 

H = H f a ree + H[ ree + H f / ee + H f B ree + H%X + H$, (8) 

does not comply with Eq. Q . That this is so follows im¬ 
mediately from the interaction terms H" 4 + , which 

under Pa b ®^E transform into + H'" B , thus prevent¬ 
ing the invariance of H under the exchange a <->• b. 

These results show that an open bipartite system pre¬ 
serves the symmetry under the exchange of its (nonin¬ 
teracting) parts a and b if and only if there is a common 
environment so that the evolution is global (nonlocal). 
This is specially relevant when studying decoherence pro¬ 
cesses in identical particle systems, a matter that will be 
analysed below, in relation with systems of two identi¬ 
cal fermions. In fact, in the particular case when a and 
b are indistinguishable subsystems, clearly the Hamilto¬ 
nian 0 can be ruled out from the start, since there is no 
possibility of distinguishing them through an interaction. 

III. PRELIMINARIES 

A. Entanglement in systems of identical fermions 

Consider that a and b represent two identical fermions. 
(Though indistinguishable, in what follows we will of¬ 
ten use the notation a and b to refer to “one” and “the 
other” fermion. Such notation is introduced for clarity 
purposes, and must not be understood as a labeling that 
distinguishes between the two fermions). Let \<p K ) and 
\4> k ') be two single-fermion states. The antisymmetric 
combination 

| V’fe ) = ^(1 I Mb - I </>*') a I Mb)- ( 9 ) 
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defines what is called a Slater determinant (and is said 
to have Slater rank 1). A composite system consisting 
of two identical fermions is regarded as separable (i.e., 
non-entangled) if and only if its density matrix is of the 
form [TT] 


PIT = J2P k I Vf) (Vf | » 


( 10 ) 


with = 1- That is, a pure state of two identical 

fermions is simply a single Slater determinant, whereas 
mixed separable states are those that can be expressed 
as a statistical mixture of pure states of Slater rank 1. 
Here, by “entanglement” we mean entanglement between 
particles (as opposed to entanglement between modes). 
Comparison of Eq. © with Eq. ( [To] ) indicates that in or¬ 
der to describe non-separable states of indistinguishable 
fermions we need to resort to more general basis {| '/T ) s } 
that include elements different from Slater determinants. 

Since there are n(n— l)/2 different k's and n = 2s+ 1, 
with s being the spin of the particle, the dimension of H _ 
equals s(2s + l). For s = 1/2 the basis {\ip ^) s } possess a 
single element, it thus possess Slater rank 1, and hence no 
entanglement is present. Therefore the fermion system 
of lowest dimensionality exhibiting the phenomenon of 
entanglement corresponds to s ^ 3/2, for which n ^ 4 
and dim"H_ > 6. Denoting with {]*)} = {11), |2) ,..., |n)} 
an orthonormal basis of the n-dimensional Hilbert space 
of each subsystem, we can identify each |i) with the states 
|s, ni s ) , with m s = — s ..., s, [7] so that 


{I 1 ) = |s,s) i |2> = |s,s- 1) 


*) = |s,-s}}. ( 11 ) 


Within this angular momentum representation, the an¬ 
tisymmetric joint eigenstates {|j, m), —j<m< j, 0 < 
j < 2s} of the total angular momentum operators J z and 
J 2 constitute a natural basis {j^T) ? } for the Hilbert 
space associated with the pair of identical fermions. 
Such antisymmetric states are those characterized by an 
even value of the quantum number j |25] 2T> . In what 
follows the notation | j, m ) is always meant to refer to 
the angular momentum representation. 


The following is a list of the antisymmetric total angu¬ 
lar momentum eigenstates for two fermions of spin | 
with the value for the concurrence (see equation ©) 
indicated on the right (we use a compact notation ac¬ 
cording to which, for instance, the ket | 0 , 0 ) stands for 

I j = 0 ,m = 0 )): 



C 

IO = | 2 , 2 ) = ^| 12 )-^| 21 ) 

0 

l^-) = | 2 ,l> = ^|13>-^|31> 

0 

|VT> = 12, 0) = 5 23) + § 14) - i|41) - 5 32) 

1 

|C) = |2,-l) = ^|24)-i|42) 

0 

k"> = | 2 ,- 2 ) = ^|34)-^|43) 

0 

1^6 > = 10,0) = | 32) - g 23) + 2 14) - 2 141) 

1 


Notice that the states |0,0) and 12,0) are maximally 
entangled, while all the other states in the list cor¬ 
respond to single Slater determinants thus have zero 
entanglement. 

Necessary and sufficient separability criteria for pure 
states of two identical fermions have been formulated in 
terms of appropriate entropic measures evaluated on the 
single-particle reduced density matrix (see [18] and ref¬ 
erences therein). For mixed states, however, the devel¬ 
opment of entanglement criteria and entanglement mea¬ 
sures remains largely unexplored. Only for fermionic sys¬ 
tems described by a single-particle Hilbert space of di¬ 
mension 4 a closed analytical expression for the amount 
of entanglement, or concurrence C(p a b), in a general 
(pure or mixed) two-fermion state p a b is known [7], 


C(pab) = max{ 0 , Ai — A2 — A3 — A4 — A5 — Xq}, ( 12 ) 

where the A/s are, in decreasing order, the square roots 
of the eigenvalues of p a bpab with p a b = Op^O -1 , where 
D is given by 


/0 0 0 0 1 0 \ 

0 0 0-100 
0 0 1 0 0 0 

0-10 0 00 
1 0 0 0 0 0 

\0 0 0 0 0 1 / 


(13) 


k is the complex conjugation operator and D is expressed 
with respect to the total angular momentum basis in the 
following order | 2 , 2 ), | 2 , 1 ), | 2 , 0 ), | 2 , — 1 ), | 2 , — 2 ), and 
*| 0 , 0 ). 


B. Entanglement between one fermion and the rest 
of the system 

Considering that the tripartite system consisting of two 
identical fermions and the environment £ is in a pure 
state \ij}), we now look for a quantitative indicator of the 
amount of entanglement between one of the fermions and 
the rest of the system. In order to do so we first notice 
the following: 

• Since both fermions are identical, the amount of 
the entanglement exhibited by one of them with 
the rest of the system must be the same for both 
of them. 

• If each fermion can be regarded as disentangled 
from the rest, then the fermions pair as a whole, 
is disentangled from the environment E. 

It is natural to use as an indicator of the amount of 
entanglement of one fermion with the rest an entropic 
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measure evaluated on the single-fermion reduced density 
matrix pf (here / is either a or b ), obtained after tracing 
the full state p = \%!)){ip\ over one fermion and over the 
environment. Let us consider the von Neumann entropy 


SvN(Pf) = -Tr(pf\np f ). (14) 


Now, the global state | if>) can be expressed in terms of 
the Schmidt decomposition 


w = ( 15 ) 

i 

where \ are the Schmidt coefficients such that )Tb = 
1 , and {\4>Y)} and {|£)}} are sets of orthonormal states 
belonging, respectively, to and He- Let pfi denote 
the single-particle density matrix corresponding to the 
two-fermion state | </>“). Then, 


Pf = J2\iPfi. (16) 

i 

Using the concavity of S v n we get 


S V N{Pf) > T, XjS vN (pfj). 


(17) 


On the other hand, it holds that m 
SvN(.Pfi') ^ In 2. 

Therefore, combining and ( fl8] ) we obtain, 
SvN^Pf) ^ yx^Mp^ a in 2. 


(18) 


(19) 


The equality sign in (18) holds only if each \(j>~) corre¬ 
sponds to one-single Slater determinant. The equality in 
( p~7] ) occurs only if all the pfi are equal to each other [27]. 
Combining these two conditions it follows that the lower 
bound in (19), S V N(pf) = In2, happens only if the two 
fermions are disentangled from the environment and dis- 
entagled from each other. That is, when the two-fermion 
state is described by a Slater determinant. 

The above considerations show that the lower bound in 


(19) corresponds to the physical situation in which each 


of the fermions has to be regarded as disentangled form 
the rest of the system. Also, the quantity 


£ = S vN (p f ) - In2, 


( 20 ) 


provides a useful quantitative indicator of the amount of 
entanglement between one fermion and the rest of the 
system. Notice that e vanishes if and only if the two 
fermions are disentangled from the environment and dis- 
entagled from each other. The measure e has a non-zero 
value if the two fermions are entangled with each other, 
or entangled with the environment, or both. 


The measure e is different from a measure of entangle¬ 
ment between the two fermions, and is also different from 
a measure of the entanglement between two fermions (as 
a whole) and the environment. In addition, e is fully 
consistent with previous approaches to entanglement be¬ 
tween identical fermions (see, for instance mi). In par¬ 
ticular, it takes into account the fact that the minimum 
correlations required by anti-symmetry do not contribute 
to entanglement. 

Basically the same reasoning used above can be applied 
to argue that 


\ ~ Mp)) ( 21 ) 

is an appropriate quantifier of the entanglement between 
one fermion and the rest of the system. This is based on 
the fact that 1 —Tr(pj) is a concave functional of pf, and 
that for a pure state | (j>) of two fermions we always have 
1 — Tr (p'j) > \ with equality if and only if |<\>) is a Slater 
determinant (16) . Consequently, we can resort to 

c -'“= yyyy ^ 

to quantify the entanglement between one of the fermions 
and the rest, when the two-fermion plus environment sys¬ 
tem is in a pure state. Here d > 4 is the dimension of the 
single-fermion Hilbert space, and the factor 2 d/{d — 2) is 
introduced so that C a \Eb lies between 0 and 1. 


C. Entanglement between one fermion and the 
environment 

For the tripartite system abE , let us consider the ob¬ 
servables of the form, 


0= ^y(Ai® I + I® Ai') (g) Bi. (23) 

i 

In this equation the A^s act on the single-fermion Hilbert 
space, that is, they correspond to observables represent¬ 
ing properties of one single fermion. I is the identity 
operator acting on the single-fermion space, and the TVs 
are observables referred to the environment. The expec¬ 
tation value of O in the tripartite state p reads 


(O) = Tr 


= Tr 


^ (a* (g) I + I (g> A;) <g> Bi 

i 

y_ Aj (8) Bi j pfE 


(24) 


where pfE is the density matrix obtained after tracing 
the global density matrix over the degrees of freedom of 
one fermion, e.g., p a E = Tr^p. 
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The observables of the form (23) are those representing 
properties referred to one fermion and the environment. 
Equation (24) means that, as far as these observables 
are concerned, all the statistics associated with quantum 
measurements are described by the reduced density ma¬ 
trix pfE■ That is, pfE jointly describes single-fermion 
and environment features of the system, including the 
concomitant correlations (both quantum and classical) 
between single-fermion properties and environment prop¬ 
erties. Consequently, it is physically meaningful to re¬ 
gard the entanglement of the state pfE (measured in the 
usual sense when considering distinguishable subsystems) 
as describing the entanglement between one fermion and 
the environment. Therefore, the entanglement between 
one fermion and the environment can be operationally de¬ 
fined as the effective entanglement between the fermions 
and the environment when only single-fermion properties 
can be measured. 

As we are dealing with pure global (tripartite) states p , 
the reduced density matrix pf E will in general be a mixed 
state of an n-level and an m-level system, m being the di¬ 
mension of He, whence we will use the negativity [251125] 
as an indicator of entanglement between one fermion and 
the environment. The negativity A f is given by the sum 
of the negative eigenvalues at of the partial transpose 
(with respect to either E or /) of the matrix PfE, 


indeed an entanglement measure whenever the total state 
is pure, i.e., described by a vector \ip) a bE [30] • As for the 
entanglement between one fermion and the environment, 
we will proceed as explained in Section |IIfC| and use 
the negativity (25) to detect the entanglement between 
/ = a, b and E. 

In the study of decoherence processes, the Kraus rep¬ 
resentation introduced in Section[H]is particularly useful, 
since it allows to represent the unitary evolution of the 
fermion (S') plus environment (E) system by the quan¬ 
tum map pQ 






Two paradigmatic quantum channels, widely used in 
studying the decoherence in open qubit systems, are 
the Amplitude Damping Channel (ADC) and the Phase 
Damping Channel (PDC). The former represents the dis¬ 
sipative interaction between the qubit and its environ¬ 
ment, and the later can represent the coupling of the 
system to a noisy environment [2]. Here we will gener¬ 
alize the main features of these channels to extend the 
corresponding quantum map to the 6 -dimensional joint 
Hilbert space of the two fermions. 


A. Amplitud Damping Channel 


M = 'Z,\a i \. (25) 

i 

By virtue of the PPT criterium, a positive value of Af 
indicates that the state pfE is entangled. 


IV. DECOHERENCE PROCESS IN 
TWO-FERMION SYSTEMS 


In this Section we shall analyse the dynamics of en¬ 
tanglement in an open system consisting of two identical 
fermions. According to the discussion in Section [TTJ we 
focus on fermions that share a common environment, and 
that evolve under different decoherence processes when 
the initial state has the form 0 . We will compute ana¬ 
lytically the entanglement between different parts of the 
complete (fernrion+environment) system, restricting the 
study to fermions with a single-particle Hilbert space of 
dimension four, immersed in a two-level environment E 
whose states are | 0 ), 11 ). 


In particular, we will resort to Eqs. (12) and (22) to 
compute the entanglement between the fermions [C a b), 
and the entanglement between one fermion and the rest 
of the system (C a \ E b), respectively. In addition, we will 
use the expression 


C E \ ab = ^2(1-Trp|) = ^2(1-Tr fa), 


(26) 


to calculate the entanglement between E and the 
fermionic subsystem. Recall that equation (26) provides 


The main feature of the ADC is that it preserves the 
total (system plus environment) excitation number. Con¬ 
sidering that in the bipartite states | 2 ,m), m stands for 
an excitation that can be exchanged with the environ¬ 
ment, the AD map in this case reads 

|2,to) s |0) b -4 sjl — p\2,m)s\0) E + y/p\2,m- 1)s|1)e 

m = — 1 , 2 , 

|0,0) S |0) B -a |0,0)s|0) B , 

|2,-2) s |0)| b -> |2, —2}s|0)_e, (28) 

where p £ [ 0 , 1 ] is a continuous parameter characterizing 
the evolution. Let us consider the initial state 

\iK0))se = \2,0)s\0)e' (29) 

According to Eq. ( [28] ) , the whole tripartite system 
evolves to 

\^{p))se = \/l -p|2, 0 )s|0)b + v^|2, —1 )s|1)_e- (30) 

We obtain the following expressions for the squared con¬ 
currences (tangles) as a function of p: 

CUP) = (1 ~P)\ (31) 

C E\ab(p) = 4 p(!-p) (32) 

Ca\Eb(p) = 1 ~ P 2 , (33) 

and plot them in Fig. [T] The solid (orange) line shows 
a typical feature of decoherence processes: as a result of 
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the interaction of the fermionic system with the environ¬ 
ment, C 2 b decreases nronotonically until its completely 
disappearance. The dotted (purple) curve, representing 
Cg\ ab , shows that along the evolution the environment 
gets entangled with the fermionic system, disentangling 
from it only at p = 1. It is also observed that the bi¬ 
partite entanglement between a and the rest ( Eb ) (green 
dashed curve) decreases at a slower rate than Cf lh , with 
C®\ Eb ^ C 2 b . I n the inset of Fig. fll we plot the evolu¬ 
tion of the negativity M{p{e)- Sucn quantity is positive 
for all p £ ( 0 , 1 ), indicating that in this interval there 
exist entanglement between each single fermion and the 
environment. 


The dynamics of entanglement induced by the ADC 
has been previously studied in the context of two initially 
entangled (distinguishable) qubits, q\ and < 72 , when only 
<72 interacts (locally) with its environment E 2 (see, e.g., 
Ref. [3T]). Though in the present fermionic system the 
environment is common to both particles so the evolution 
is non-local, the comparison between the distinguishable- 
qubit and the identical-fermion case seems useful to 
evince the main features that distinguishes one and the 
other type of evolution. In particular, in the qubit case, 
it is found inn that as a result of the decoherence chan¬ 
nel, there is an entanglement swapping between qiq 2 and 
q\E 2 , that is, the initial (p = 0) entanglement between 
<71 and q 2 is completely converted (at p = 1 ) into entan¬ 
glement between <71 and E 2 . In this sense the net effect 
of the ADC is to redistribute and transfer the initial bi¬ 
partite entanglement without looses. In the fermion case 
this no longer holds. This can be seen by taking p = 1 
in the state (301 and observing that since 12, — 1) is a 
two-fermion separable state (see the Table), the final tri¬ 
partite state |V , ( 1 ))s_b = | 2 , — 1 )s| 1 )b is completely dis¬ 
entangled. Thus the initial (maximal) entanglement be¬ 
tween the fermions is finally lost due to the decoherence 
process, yet during the evolution (i.e., for 0 < p < 1 ) the 
entanglement redistribution due to the ADC is of course 
present, as seen in Fig. [T] These observations would 
thus indicate that the open qubit system is more robust 
against decoherence than the identical-fermion system. 
This is reinforced by the fact that for an initial max¬ 
imally entangled state in the qubit case we have {31, 
Cg iq2 (p) = 1 —p, whereas here the tangle between the two 
fermions is given by Eq. (31), namely C^ b (p) = (1 — p) 2 , 
so the entanglement between the fermions decays faster 
than the entanglement between the qubits. 

As a second example we consider the initial state: 


|^(0))sb = (a|2, l)g + /3|2, -l)s)|0) E , (34) 


with |a | 2 + |/?| 2 
evolves to 


1. Applying the map (28) the state 


I 4>{p))se = a[\/l - p|2, l)s|0) B + y/p\2, 0) s |0)_e] (35) 
+ /?[\A - Pl 2 , -l)s|0)s + Vp|2, -2)s|1).e], 



FIG. 1: (Color online). Entanglement evolution for the initial 
state ( |29| ), under Amplitud Damping decoherence channel. 
Evolution of the entanglement between the fermions (orange 
solid line), evolution of the entanglement between the envi¬ 
ronment and the pair of fermions (violet dotted line), and 
evolution of the entanglement between one fermion and the 
rest of the system (green dashed line). Inset: evolution of 
the negativity of the reduced state of one fermion and the 
environment. All depicted quantities are dimensionless. 


and in this case we obtain: 

Cabb) = [(1 -p)(aP* +a*!3) -p|a| 2 ] 2 , 

C|| abiP) = 4 P( 1 -P), 

Cl lEb ( P ) = 2 — 4p|a| 2 |/3| 2 — [1 + (1 — p) 2 ]|a| 4 

- 2[p + (1 — p) 2 ]|/3| 4 . (36) 

Three qualitatively different cases will be now analysed 
for different values of the parameters a and /?. For a = 0 
(/3 = 1 ) we get C% b (p) = 0 , thus the fermions remain 
in a separable state along the whole evolution, whose 
only effect is that of continuously transforming the state 
12 , — 1 ) (at p = 0 ) into | 2 , — 2 ) (at p = 1 ) without modify¬ 
ing the entanglement between the identical parties. For 
a = /3= 1/V2, we obtain the following: 

Cab(p) = (i - , 

c 2 E \ab(p ) = 4p(i-p). 

Ca\Eb(p) = 1 — ^ P 2 ‘ ( 37 ) 

Thus C 2 h (p) decreases monotonically from its maximum 
value to zero, at p = 2/3, where it starts to increase 
as shown in Fig. [2] Such increase in the entanglement 
between noninteracting entities is a result of the non¬ 
local dynamics due to the common (or global) environ¬ 
ment mm- A more drastic example of the increase in 
the entanglement between the fermions due to the global 
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FIG. 2: (Color online). Entanglement evolution for the ini¬ 
tial state (341 with a = (3 = 1 /\/2 under Amplitud Damping 
decoherence channel. Evolution of the entanglement between 
the fermions (orange solid line), evolution of the entanglement 
between the environment and the two fermions (violet dotted 
line), and evolution of the entanglement between one fermion 
and the rest of the system (green dashed line). Inset: evolu¬ 
tion of the negativity of the reduced state of one fermion and 
the environment. All depicted quantities are dimensionless. 


environment can be seen by taking /3 = 0 (a = 1) in 
Eqs. (361. In this case the initial state |2,1) is sepa¬ 
rable (see Table), but as p increases the entanglement 
between fermions increases as well, since C 2 b {p) = p 2 , 
so that at p = 1 the fermions end up maximally entan¬ 
gled. In other words, when applied to appropriate initial 
states, the decoherence channel is capable of increasing 
the entanglement between the fermions. 

The negativity for the case a = /? = l/\/2, shown 
in the inset of Fig. [2j is qualitatively the same as in 
the previous case (initial state ( |29| ) subject to the ADC). 
Again, in the interval p £ (0,1) there exist a nonzero 
entanglement between one fermion and the environment. 


B. Phase Damping Channel 

This process describes the loss of quantum information 
with probability p without any exchange of energy. The 
PDC is described by the quantum map 


\j,m) s \0) E —t y/l-p\j,m)s\0) E + Vp\j,m)s\V)E 
j = 2,?n = -2,..., 2 

|0,0)s|0)_e -> |0,0) S |0) B . (38) 

We apply the map to the initial state 

\v(0))se = (<5|2, 0) s + *7|0, 0)s)\0)e, (39) 

with |<51 2 + |y| 2 = 1. The evolved state reads 

\V(P)) = 5 |2,0) s | P(p)) E + *7 |0, 0) s |0) B , (40) 


where we defined 


\p(p))e — P I®)b + Vp |1)e • (41) 

The squared concurrences for this case are given by 
Cab(p) = C(p,5, 7 ) 

- VC 2 (Pi 7) - [C(P^,7) - 2p|^| 2 |7| 2 ], 
C E\ab(p) = 2[1 — C{p, <5, 7 )], 

Cl\ Eb {p) = 1-(1-p)(<5*7-^7*) 2 , (42) 

where C(p,S,j) = |5| 4 + I 7 I 4 + 2(1 — p)|< 5 | 2 | 7 | 2 . Setting 
S = = — iS , we get the following expressions: 

c lb(.P) = 

C E \ab(P) = Pi 

C a \Eb(p) = 1; (43) 

and plot them in Fig. [ 3 ] Now C 2 b decreases linearly in 
p whereas C%i ab increases at the same rate, so the sum 
C 2 b + C^ ab remains constant along the evolution, and 
equal to C 2 ^ Eb . Unlike the previous (ADC) case, here 
the environment ends up being maximally entangled with 
the bipartite system S at the expense of the disentangling 
between the fermions. 



FIG. 3: Entanglement evolution for the initial state (391 with 
S = l/\/2 ,7 = — iS under Phase Damping decoherence chan¬ 
nel. Evolution of the entanglement between the fermions (or¬ 
ange solid line), evolution of the entanglement between the 
environment and the two fermions (violet dotted line), and 
evolution of the entanglement between one fermion and the 
rest of the system (green dashed line). All depicted quantities 
are dimensionless. 


As for the entanglement between one fermion and the 
environment, we resort to Eq. ( |40| ) to obtain the reduced 
density matrix p aE (p) = T r b \p(p)) (r?(p)|. Direct calcu- 





























lation leads to 


Pa E ( P ) = ^[|1) (l| a + 14) (4| a ] \(fi + (p)) (lfi + (p)\ E 

+ ^[|2) (2| a + |3) (3| a ] (ip-(p)\ E ,(44) 

where \i) a (i = 1 ,2,3 ,4) are the single-fermion states 
defined in Section III A[ and \ip± (p)) stands for the vector 


I <P±(p)) = ^[|0)±|P(p))] 


= ^[(1± V / 1^)|0)±Vp|1)]. (45) 


Equation (44) shows that the state p a E is separable for all 


p, so that any measure C a E quantifying the entanglement 
between one fermion and the environment vanishes, i.e., 
Co.e{p) = 0. 

It is interesting to observe that this latter result, to¬ 
gether with Eqs. (43), coincide with the qubit concur¬ 


rences obtained in the previously discussed 3-qubit sys¬ 
tem m- Thus, contrary to what happened in the ADC 
case, the effect of the PDC on both (qubit and fermion) 
systems seems to be the same regardless of the local (or 
non-local) nature of the interaction. In order to go fur¬ 
ther in the comparison between the distinguishable-qubit 
and the identical-fermion case, we recall the monogamy 
inequality 


C*\ jk - Cl -C 2 k > 0 (46) 

satisfied by the usual concurrence, i.e., involving distin¬ 
guishable qubits i, j and k [321 - Motivated by this in¬ 
equality we define 

R a = C 2 a]Eb - C 2 aE - C 2 ab (47) 


genuine tripartite entanglement of those states pertaining 
to the GHZ-type family j34j. Thus, we can say that r(p) 
here measures the genuine tripartite entanglement shared 
by the two fermions and the environment, and that the 
state (401 is the analogous of the GHZ-family states for 


systems involving two fermions and their common envi¬ 


ronment. According to Eq. (49), the genuine tripartite 
entanglement increases linearly in p , and becomes max¬ 
imum at p = 1, where the ‘fermion-environment GHZ’ 
state, 

I GHZ) se = 1*7(1)) = -^(|2, 0) s |1) B + |0,0) s |0) B ), 

(50) 

is reached. The state (50) shares with the usual (3-qubit) 


GHZ state the property of having maximal genuine en¬ 
tanglement, while having zero entanglement between the 
parties when one of them (any) is traced out. 


V. ENTROPIC ENTANGLEMENT CRITERIA 

As we mentioned before, the particular case of systems 
of two identical fermions with a four-dimensional single¬ 
particle Hilbert space (the simplest fermion system ad¬ 
mitting entanglement) is the only one for which we have 
a closed, analytical expression for the concurrence. No 
such expression is known for fermion systems of higher 
dimensionality. 

In order to study the entanglement dynamics of sys¬ 
tems of ./V-fermions undergoing decoherence it is possible 
to use an entanglement indicator based upon entropic cri¬ 
teria m- In this section we will use our previous results 
for the case of systems of dimension 4 to investigate the 
efficiency of these criteria. 

All separable states (pure or mixed) of N identical 
fermions comply with the entropic inequalities 


and 


S { *\p F ) + lnN>s£\p f ), (51) 


Re = C% ab - C 2 Ea - C 


Eb'> 


(48) 


where C a E = Ceo. = C-Eb is an appropriate measure (con¬ 
sistent with the previously defined concurrences) of the 
entanglement between the two-level environment and the 
four-level fermion. Using C{e(p ) =0 and Eqs. (431 leads 
to 


R a (p) = Re(p) = r{p) =p^0. (49) 


where S^ is the Renyi entropy with a > 1, pf is the 
global ./V-fermions density matrix, and pf is the single¬ 
particle reduced density matrix. The equality sign in the 
above inequality occurs, for instance, in the case of pure 
separable states. Now, if ones considers the quantity 

Q ia) (pf) = (pi) - 4 a) (pf) - In N, (52) 
then for all separable states (pure or mixed) one has 


Since R a ,E encodes information of the entanglement 
that cannot be written as entanglement between two par¬ 
ties (hence reflect multipartite entanglement), a positive 
value of R a ,E exhibits the presence of tripartite entan¬ 
glement. Moreover, since R a = Re-, such tripartite en¬ 
tanglement is the same in all bipartitions (fermion(rest, 
environment|rest). In the 3-(distinguishable) qubit sys¬ 
tem the corresponding residual entanglement is just the 
3-tangle T t j k = r = C 2 ,- k - Cfj - C 2 k , which measures the 


Q (a \p F )< 0 . (53) 

Therefore for Q > 0 one knows for sure that the state 
is entangled: 

Q ( ' <x \pf) > 0 => pf entangled. (54) 

Since the converse (p E entangled => Q( a > > 0) does not 
hold in general, the condition Q > 0 detects some 
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(mixed) entangled states but not all of them. For two- 
fermion states the entanglement criterion improves as a 
increases and is the most efficient in the limit a —» oo El- 
Note that the criterion associated with the von Neumann 
entropy constitutes a special instance, corresponding to 
the particular value a —> 1 of the Renyi entropic param¬ 
eter. 



FIG. 4: Entropic entanglement indicator evolution (dashed 
line), and evolution of the fermionic concurrence (solid line) 
for the initial state (291 under ADC. Inset: Evolution of the 
entropic entanglement indicator (dashed line) and fermionic 
concurrence dynamics (solid line) for the initial state (341 
with a = /? = l/%/2 under ADC. The logarithms in the en¬ 
tropic entanglement indicator are taken to the base 2 and all 
depicted quantities are dimensionless. 


The quantity Q^ can thus be regarded as an entan¬ 
glement indicator whose evolution under decoherence can 
be investigated. In particular, we study the evolution 
of Q under decoherence for two previously considered 
examples setting a = oo. In this case we have (with 

Q (oo) = Q) 

Q(Pf) = (Pf) - <S'k° ) (, Pf ) - In2, (55) 


with 


s r°\pf) = -lnA„ 


(56) 


where A max is the largest eigenvalue of pp- Figure [4] 
shows the evolution of Q (dashed lines) for the initial 
state (29) under the Amplitude Damping decoherence 


channel. It is observed that the evolution of the entan¬ 
glement indicator is qualitatively the same as the evolu¬ 
tion of the fermionic concurrence (solid lines). The re¬ 
semblance is stronger in the case of the initial state (34) 


with a = /3 = 1 /x/2 ( inset). However, in any case, we 
can conclude that Q is a reasonably good entanglement 
detector. 

The entanglement indicators Q^ a \ as well as the en¬ 
tanglement measures considered in Section III, are not 


straightforwardly measurable, in the sense that they are 
not equal to (or function of a small number of) expecta¬ 
tion values of quantum mechanical observables. However, 
if the global state of the two-fermion system under con¬ 
sideration is first reconstructed via appropriate quantum 
state tomography techniques, then the aforementioned 
entanglement quantities can also be experimentally de¬ 
termined. More directly measurable entanglement indi¬ 
cators for fermionic systems have not yet been as inten¬ 
sively investigated as those for systems with distinguish¬ 
able parts. However, some progress in this direction has 
been done. For instance, the entanglement indicators 
advanced in |20| . based upon uncertainty relations, are 
expressed in terms of expectation values of measurable 
quantum observables. 


VI. CONCLUSIONS 

We studied the effects arising from the interaction be¬ 
tween a quantum system of two identical fermions and 
the environment. We showed that for the exchange sym¬ 
metry to be preserved, the evolution of the system must 
be global, or nonlocal, in the sense that each fermion in¬ 
teracts separately with a common environment. Thus, in 
order to analyse the dynamics of entanglement under two 
paradigmatic decoherence channels widely studied in the 
context of local qubit dynamics, we generalised and ex¬ 
tended the Amplitude Damping Channel and the Phase 
Damping Channel to the joint Hilbert space of the two 
fermions. 

In order to achieve a more complete analysis of the 
evolution of the entanglement in the tripartite system 
(fermion+fermion+environment), it was necessary to de¬ 
fine two measures of entanglement: one that quantifies 
the entanglement between one fermion and the rest of the 
system (fermion+environment), and one that quantifies 
the entanglement between one (any) of the fermions and 
the environment. With these tools we were able to study 
the dynamics of entanglement for some initially entan¬ 
gled states subject to the ADC and the PDC. Compar¬ 
ison with the 3-qubit case was made, and new insights 
into the mechanism of entanglement evolution in open 
systems of identical particles were revealed. 

In the case of the Phase Damping Channel, and by re¬ 
source of a monogamy relation, we were able to detect 
genuine tripartite entanglement and determined an anal¬ 
ogous to the GHZ-state involving the two fermions and 
the environment. Further progress in relation with tri¬ 
partite entanglement and monogamy relations in these 
kind of systems is, however, inherently constrained by 
the advance in the problem of quantifying entanglement 
in systems involving identical particles, and multipar¬ 
tite systems in general, a problem that remains far from 
closed. 

Finally we showed the dynamics of an entanglement in¬ 
dicator based on an entropic criteria which can be used 
to study decoherence in more general (higher dimensions) 
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systems of identical fermions. For the four-dimensional 
case studied here, the entanglement indicator turned out 
to be a reasonably good indicator of entanglement be¬ 
tween the pair of fermions. 

A possible experimental scenario in which to consider 
the kind of processes discussed in the present work could 
be provided by a system consisting of two electrons 
in two laterally coupled quantum dots [Bj. This sys¬ 
tem allows for the implementation of quantum informa¬ 
tion related processes such as quantum gates, and can 
be described in terms of an effective four dimensional 
single-particle Hilbert space (leading to a six-dimensional 
two-fermion Hilbert space). The relevant single-particle 
Hilbert space is spanned by single-electron states that, 
in self-explanatory notation, can be expressed as {|A 'f' 
),|A 4-), \B I)}, where {|A),|H)} denote two or¬ 


thogonal electronic spatial wave functions (orbitals) pre¬ 
dominantly located around two particular locations in 
the double quantum dot. This two-electron double quan¬ 
tum dot system then constitutes a possible experimental 
realization (of technological significance) of the type of 
two-fermion systems considered in the present contribu¬ 
tion. 
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